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We calculate axial vector and vector form factors for the semileptonic decay Dj — ^ using HISQ 
valence quarks on MILC ensembles with 2 + 1 flavours of asqtad sea quarks. Using twisted bound- 
ary conditions to tune the quarks' momenta, we compute form factors at = on coarse lattices. 
We find V (0) = 0.903 (67), A i (0) = 0.603 (20), A2 (0) = 0.401 (80) and Aq (0) = 0.686 (17), which 
we compare to experimental data and previous quenched lattice QCD calculations. 
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1. Introduction 

The decay rate for the semileptonic meson decay -> 0/v contains nonperturbative QCD 
form factors, about which we obtain information from lattice QCD. The quark level transition is 
c ^ s and the decay is from a pseudoscalar meson to a vector. 

Labelling the pseudoscalar as Dg with momentum p and the vector meson as with momentum 
p' and polarisation vector e, the general form for a pseudoscalar to vector matrix element is [|l[]: 

W,e)\V^ -A^\D,{p)) = ^:!ft^e'p-/v{q') (1.1) 

- {mD,+m^)e^Ai{q^) 

+ ^^{p+prMq') 

+ 2m^^q^A3iq^) 

r 

-2m0— r-(7'^Ao(<?^) 

where we have the relations 

A3(.^) = ^^^A,(,^)-!!!^A.(,^) (1.2) 

and 

A3(0)=Ao(0). (1.3) 

V{q^) is the vector form factor. There are three independent axial vector form factors as AT,{q^) is 
just a combination of Ai{q^) and A2{q^) and at = there are only two. Calculating the form 
factors at = allows us to compare to experiment for the experimentally accessible A\{q^), 

Am , _ v(o) 



A2{q^) and V{q^) form factors. BaBar quote results for the quantities r2 = and ry — 



2. Methods 
2.1 Lattice setup 

To compute semileptonic form factors, we calculate strange quark propagators using a random 
wall source on a timeslice at x. A strange quark propagates to the sink, y, and then is used to 
make an extended propagator in which a charm quark propagates back to the current insertion at 
z, which is summed over a timeslice. At the current insertion, the extended propagator is joined to 
the strange propagator to make the 3 point function. 

We use HISQ [Q] valence quarks, which means that suitable phases for staggered quarks have 
to be inserted at the source, sink and current insertion. The gauge configurations used are from 
MILC ensembles [^, which include 2+1 flavours of asqtad sea quarks. Currently we have per- 
formed calculations on only one ensemble, which has j3 = 6.76 and V = 20^ x 64, so a ~ 0.12fm. 
In general, the 3 point correlation function is 

{W)yiiW) W)y5yvW) viy)r5viy')), (2-i) 

where we allow for point splitting at the source, sink and current insertion. 
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(a) Dispersion relation. 



(b) Ground state amplitude. 



Figure 1: To check that introducing twisted boundary conditions gives us a physical momentum, we can 
calculate the speed of light for various mesons and momenta, = ^ (Fig. 1(a) ) and look at the ground 



state's amplitude. A, as a function of meson energy. We expect A v£ to be constant (Fig. 1(b)). 





l-link 


local 


experiment 




1086(4) MeV 


1083(4) MeV 


1019 MeV 


U 


248(4) MeV 


233(13) MeV 


227(4) MeV 




Z= 1.095 


Z = 0.97 


PDGr,+, 



Table 1: The masses and decay constants for our ^ mesons calculated using the l-link and local vector 
operators. We include experimental values from the PDG and note the lattice normalisation factors used 
for each calculation. 



2.2 Twisted boundary conditions 

To calculate the form factor at = 0, we use a with non-zero momentum and a Ds meson at 
rest. We use twisted boundary conditions where we insert a tunable phase, Q, at the boundary 
so ^|f{x +L) = \if{x)e'^ . In Figure [T], we check that this introduces a physical momentum by looking 
at the dispersion relation and the amplitude of the ground state. 

2.3 (j) 2 point correlators 

To extract matrix elements from 3 point correlators, we fit the 2 and 3 point functions simulta- 
neously. The 2 point functions give us the masses and decay constants of the ^ and mesons in 
the 3 point amplitudes. The mass is known to be accurately calculated with HISQ quarks 
We also extract and f^, using either a taste-singlet l-link or local vector operator. Table [I] gives 
the masses and decay constants of the 1-hnk and local (j) mesons, which are in agreement. 

The meson is not gold-plated, which means that it can decay to KK with a small width, an 
effect we do not fully include. We have only done these calculations at one lattice spacing, so we 
cannot yet tell how much of the difference between the calculated lattice and continuum <p masses is 
caused by discretisation effects. With results on finer lattices, the effects of the non-gold-platedness 
of the can be studied in more detail. The normalisation factor for the l-link vector operator is 
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from Ref. [10| and a perturbative normalisation is used for the local vector |11]. We extract 

1, using 



from the PDG results for Fg+g 



471 



(2.2) 



2.4 Operator normalisation 



2.4.1 1-linkA 



To normalise the one link axial vector operator, we use simultaneous fits to 

(i//(x)75i/A(x) viy)r5rn¥iy+i^)) and {\j/{x)Y5vix) viy)r5viy)), 



(2.3) 



making use of the facts that the operators \jf{x)Y5Y{x) and \ff{x)y5YnV{x + fi-) have the same taste 
and the pseudoscalar current is absolutely normalised. 

The pseudoscalar and axial vector operators (in the continuum) can be related by 



/'^.(0|A^|Po) = (mi+m2)(0|P|Po). 



(2.4) 



The fit output gives the amplitude for both the pseudoscalar and the axial vector operators, which 
we obtain for several values of the meson momentum. The ratio of amplitudes is a straight line 
against with a gradient of the Z factor times the appropriate mass combinations. In Figure ^, 
we plot this for various charm and strange pseudoscalar mesons. For cs, the current we need for 
Ds — (j), we obtain the value Z = 1.065(5). 
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Figure 2: The ratio of amplitudes for the pseudoscalar and axial vector operators goes as so the 

Z factors can be extracted from the gradients of these lines. 



2.4.2 Local 

We can normalise the local At component by creating non-Goldstone pseudoscalar mesons to 
compare to Goldstone mesons 

{¥ix)Y5Yt¥{x) ¥{y)y5rt¥{y))- (2-5) 



4 



Axial vector fomi factors in (j) semileptonic decays from lattice QCD 



Gordon Donald 



With HISQ, the masses are very close (for the D^, the mass difference is 4.5 MeV) and the ampli- 
tudes, ai^G and ac for the non-Golds tone and Golds tone mesons respectively, are related by 

ZEd.Qi^g = + ms)aG, (2.6) 

where we have included the Z factor for the local axial current. For cs, we obtain Z = 1.036(4) for 
the local operator. Note that because the masses are so close, we can also use non-Goldstone 
Ds mesons in 3 point functions, which is useful when calculating the vector form factor, V{q^). 



3. Ds — )• form factors 

3.1 Ai{q^) form factor 

By setting up the polarisation and momentum such that £ • ^ = 0, we can isolate the Ai{q^) 
form factor: 

{^{p',e)\A^'\D,{p)) = {mD,+m^)e^'Ar{q^). (3.1) 

We can use either the l-link or local operator for the <p source, along with the appropriate axial 
vector current insertion. 

The determinations of the Ai{q^) form factor using both the local and 1-link operators are 
plotted in Figure ^, along with the other form factors for — )• ^. 

3.2 V{q^) form factor 

We can extract the vector form factor, V{q^), using a vector current insertion from Equation 
(1.1) provided at least one of the and <p mesons is not at rest. The 3 point function (with at 
rest) is 

lis 

{y{x)y^Y{x + v) mraViz) mrsYMy)) = ^e^'Eoy^Viq^), (3.2) 

mo, + 

where jJ., v and a are the three spatial directions. We use the non-Goldstone D., so need only 
include one link of point splitting to obtain an overall taste-singlet three point correlator for the 
pseudoscalar to vector transition with a local vector current insertion. 

We use a local vector normalisation of Z = 1.021(2) which is obtained [ |T0| ] by comparing a 
non-Goldstone Dg — )• rjs decay through a local temporal vector current to the Goldstone Ds — )• rjs 
with a (normalised) 1-link vector current. V{q^) is plotted in Figure^ 

3.3 Ao{q^) and A2(^^) form factors 

The form factor Ao(^^) can be extracted using the pseudoscalar density: 

{m, + ms){^{p',e)\P\Ds{p))=q^{^ip',e)\A^'\Ds{p)) (3.3) 

= Inline ■ qA(){q^) . 

The pseudoscalar form factor is not experimentally accessible, but is useful because it provides a 
route to obtaining A2{q^). 

We can use Equation ( |T3| ) at q^ = 0, A3(0) = Ao(0) to extract A2(0) from Ao(0) and Ai(0), 
although this does not hold at other values of q^. However we can calculate {^{p' ,e)\A^i\Ds{p)) 
with e -q^O and use our determinations of Ai (^^) and Ao(^^) to find A2(^^) with Equation (|0|), 
provided ^ is not too large. 
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Figure 3: The form factors for (j). Because of the kinematics, only Ai {q^) contributes at qlj^.,-^ (where 
both the Ds and (p mesons are at rest). All the form factors are calculated with p,j,a = 0.3 and p^a = 0.4195, 
where the latter corresponds to = 0. We also plot BaBar's values for the experimentally accessible Ai{q^), 
A2{q^) and V{q^) form factors at q^ = 0. 





HPQCD 2011 


BaBar 2008 [|] 


UKQCD 2001 [§ 


MO) 


0.594(22) 


0.607(11)(19)(18) 


0.63(2) 


LocalAi(O) 


0.603(20) 






V{0) 


0.903(67) 


1.122(85)* 


0.85(4) 


Ao(0) 


0.686(17) 




0.63(2) 


A2(0) 


0.401(80) 


0.463(61)* 


0.62(5) 


rv 


1.52(12) 


1.849(60) (95) 


1.35(7)* 


ri 


0.62(12) 


0.763(71)(65) 


0.98(8)* 



Table 2: Table of form factors at q^ = 0. We quote our determinations of A\{Q) using 1-link and local 
operators separately. Entries marked with * have been extracted from the published results. BaBar's results 
include statistical and systematic errors. The third error quoted on their value of Ai (0) is from theoretical 
input. 



3.4 Comparison with experiment and previous calculations 

Table ^ compares our results at = to BaBar's and those of previous quenched lattice QCD 
calculations. Our results for A2 (0) and r2 use our local A i(0) result, but the 1 -link result is in close 
agreement. We calculate ^2 directly from our results for Ai(0) and A()(0) for better control of the 
errors. 

In Figure ||, we plot the shape of the form factors as a function of q^. We expect the 
dependence to be relatively mild, and approximately the same, for all the axial vector form factors 
and steeper for the vector form factor. BaBar assume the shape of the form factor and extract its 
value at = from the total semileptonic rate. Our results for A 1 (0) and A2(0) are in agreement 
with BaBar and our V(0) result is a little lower. In the previous quenched lattice QCD results, no 
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difference was found between Ai (0), A2(0) and Ao(0). It is clear from Figure ^ that our results are 
different for these form factors at = 0. 

3.5 Future work 

To understand the discretisation effects better and also the importance of the non-gold-platedness 
of the mesons, we intend to repeat our calculations on MILC fine configurations, along with the 
determination of the Z factors and semileptonic matrix elements. We will also look at using lighter 
sea quark masses as, although neither the (/> nor the Ds contain light valence quarks, the non-gold- 
platedness of the may be an issue. 

The V{q^) form factor is also important for radiative decays, where the photon is coupled by 
a vector operator, such as — 77^7. We will calculate this with the same method used here for 

We can extend our axial vector form factor calculations to heavier quarks. Using HISQ for the 



max) 



heavy quark, we will extrapolate to the b quark mass [12] and calculate quantities such as Ai {q] 
for B-^D*. 
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